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Abstract. We describe the cross-correlation measurements being carried out 
on data from the LIGO Livingston Observatory and the ALLEGRO resonant bar 
detector. The LIGO data are sampled at 16384 Hz while the ALLEGRO data 
are base-banded, i.e., heterodyned at 899 Hz and then sampled at 250 Hz. We 
handle these different sampling parameters by working in the Fourier domain, and 
demonstrate the approximate equivalence of this measurement to a hypothetical 
time-domain method in which both data streams are upsampled. 
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1. Introduction 

Analysis is currently underway to search for the signature of a stochastic gravitational- 
wave background (SGWB) by correlating the signals of the 4 km interferometer (IFO) 
at the LIGO Livingston Observatory (LLO)m2] with the ALLEGRO resonant bar 
detector|j|l As described elsewhere[H], the LLO- ALLEGRO experiment is 

sensitive to a higher frequency range than the corresponding experiment using LLO 
and the IFOs at the LIGO Hanford Observatory (LHO)|7|, thanks to the relative 
proximity of the ALLEGRO and LLO sites. Additionally, the ALLEGRO detector can 
be rotated, changing the response of the experiment to a SGWB and thus providing a 
means to distinguish gravitational- wave (GW) correlations from correlated noise. [8]. 

The results of a cross-correlation measurement using data taken during LIGO's 
second science run (S2) will be reported in the near future. [2] The present work 
describes some details of the analysis procedure used, notably the handling of the 
different sampling rates of the data taken by the two detectors, and the heterodyned 
nature of the ALLEGRO data. The procedure applied here may prove useful in 
coherent measurements involving data sampled at different rates, such as data from 
the LIGO and Virgo[H3j IFOs. 
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dfL = n(x a x b -y a y") (2) 



2. Stochastic Gravitational- Wave Backgrounds 

2.1. Definitions 

A gravitational wave (GW) is described by the metric tensor perturbation h a b{r, t). A 
given GW detector, located at position rdct on the Earth, will measure a GW strain 
which is some projection of this tensor: 

h(t) = h ab (f dcU t)d ab (1) 

where d ab is the detector response tensor, which is 
1 

4 (ifo) - 2 ( 

for an interferometer with arms parallel to the unit vectors x and y and 

d (bar) = U U ( 3 ) 

for a resonant bar with long axis parallel to the unit vector u. 

A stochastic GW background (SGWB) can arise from a superposition of 
uncorrelated cosmological or astrophysical sources. Such a background, if isotropic, 
unpolarized, stationary and Gaussian, will generate a cross-correlation between the 
strains measured by two detectors. ^21 EH In terms of the continuous Fourier 
transform defined by 

p oo 

a(f) = / dta(t) exp(-»27r/[t-to]) , (4) 



(where to is an arbitrarily-chosen time origin) the expected cross-correlation is 

< W)M/')> = \&U f) p**U) 7i 2 (/) (5) 

where 

7i 2 (/) = d lab df J- JJ d 2 n fl P TT ™ b d e^f^-fo/c (6) 

is the overlap reduction function^J between the two detectors, defined in terms of 
the projector P T ' Tn '^ onto traceless symmetric tensors transverse to the unit vector n. 
Figured shows the overlap reduction functions for several detector pairs of interest. 

P gw (/) is the one-sided spectrum of the SGWB. This is the one-sided power 
spectral density (PSD) the background would generate in an interferometer with 
perpendicular arms, which can be seen from JSJ and the fact that the overlap reduction 
function of such an interferometer with itself is unity. Since the overlap reduction 
function of a resonant bar with itself is 4/3 (see f° r more details), the PSD of the 
strain measured by a bar detector due to the SGWB would be (4/3)P gw (/). 

A related measure of the spectrum is the dimensionless quantity fi gw , the GW 
energy density per unit logarithmic frequency divided by the critical energy density 
p c need to close the universe: 

u/)-% = S/ 3 a). (7) 



p c df 3-ff, 
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Figure 1. The overlap reduction function for LIGO Livingston Observatory 
(LLO) with ALLEGRO and with LIGO Hanford Observatory (LHO). The three 
LLO-ALLEGRO curves correspond to the three orientations in which ALLEGRO 
was operated during LIGO's S2 run: "XARM" (N72°E) is nearly parallel to the x- 
arm of LLO ( "aligned" ) ; "YARM" (N18° W) is nearly parallel to the y-arm of LLO 
( "anti-aligned" ) ; "NULL" (N63°W) is halfway in between these two orientations 
(a "null alignment" midway between the two LLO arms). The LLO-LHO overlap 
reduction function is shown for reference. 



2.2. Detection Method 

The standard method to search for such a background is to cross-correlate the outputs 
of two gravitational wave detectors ^21- If each detector signal si i2 (t) is assumed to 
be made up of a gravitational wave component h\, 2 {t) plus an instrumental noise 
component 

*i,a(*) = /n,a(*)+»i,a(*) (8) 
and the noise in the two detectors is approximately uncorrelated but significantly 
larger than the gravitational-wave signal, then the average cross-correlation should 
come from the stochastic GW background: 

W)Mf)) « (K(f)Mf)) = ^(/-/0^w(/)7i2(/) (9) 

while the average auto-correlation should come from the noise: 

(^,2(/)si,2(/')> « K 2 (/)ni,2(/0> = \s{f /') P ll2 (/) (10) 
If we construct a cross-correlation statistic 

Y c = f dt 1 dt 2 s 1 (t 1 )Q(t 1 -t 2 )s 2 (t 2 )= f dfs{U)QU)Mf) (11) 
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the expected statistics of ff 1 II) are given bv[T^| 

T 
~2 



(Y c ) 



dfl(\f\)P m (f)Q(f) 



and 



((Y c W?) « - 



dfP 1 (f)P 2 (f) Q(f) 



(12) 



(13) 



Using 1)12)1 and (|13fl . the optimal choice for the filter Q(f), given a predicted shape 
for the spectrum P gw (/) can be shown [TB] to be 

7(1/1) Pprif) 



Q(/)« 



(14) 



This has negligible support when P±(f) or P 2 {f) is large, which allows one to limit 
the integration in l|ll|l to a finite range of frequencies. 

In Sec. |3 we consider how the approximate continuous-time expressions 1)120 and 
<|13fl manifest themselves given the different discrete sampling parameters of the LIGO 
and ALLEGRO detectors. 



3. Data Analysis Technique 

3.1. Frequency- Domain Method 

As described in the frequency range of sensitivity correlation measurements 
involving the ALLEGRO resonant bar detector and the LIGO Livingston Observatory 
40 km away is determined by the sensitive bandwidth of ALLEGRO E] , which is 
limited to a band a few tens of hertz wide near 900 Hz. Due to the relatively narrow- 
band nature of the detector response, the data are base-banded, i.e., heterodyned at 
899 Hz (during the S2 run) and downsampled to 250 Hz, so that the data represent 
a frequency range from (899 - 125) Hz = 774 Hz to (899 + 125) Hz = 1024 Hz. The 
LIGO data are sampled at a frequency of 16384 Hz and therefore represent frequencies 
ranging from -8192 Hz to 8192 Hz. 

Previous work 6 proposed resampling and heterodyning the LIGO data before 
cross-correlation with ALLEGRO data. However, since this would involve not only 
downsampling by powers of two but also upsampling by a factor of 5 3 , the current 
approach is to approximate 1)11)1 in the frequency domain. 

First, we describe the relationship between the discretely sampled data in each 
instrument and the underlying continuous time series sx^it)- The first time series, 
which we will take to be the LIGO data, is sampled at a frequency {8t\)~ l , which for 
the sake of this example we will take to be 2048 Hz.J Before sampling, it is low-pass 
filtered to avoid aliasing of higher-frequency data into the analysis band. We idealize 
the effects of this anti-aliasing filter by defining 

5 (».(/) I/Ki 

1° l/l>ife 
and writing the discretely-sampled signal as 

Si\j] = S 1 (t +j6t 1 ) i = 0,...,JVi-l (16) 

| The actual sampling rate is 16384 Hz, but we digitally downsample the data before analysis, a 
process which is straightforward enough that we don't need to explicitly address it here. 
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The ALLEGRO data are first heterodyned at f% = 899 Hz to produce a continuous 
time series 

«§(t) = e~ l2 ' Tf ^ t - t "\s 2 (t) (17) 
then anti-alias filtered 

s UU) -w/f + /) i/i<A 

1° l/l > ife 

and sampled at (f^) -1 = 250 Hz to obtain 

S%[k] = S^{t + k 5t 2 ) k = 0, . . . ,N 2 - 1 (19) 

Note that this signal is now intrinsically complex. The "unheterodyned" time series 
Sift) — e l27r ^ 2 ^~ to 'S 2 it) is a band-passed version of the original signal 



r s- 2 (/) 


\f- 




^ 2 5*2 


1° 


\f- 




> 25tl 



= , ,„ 2 (20) 



While S2(*) is real, S 2 it) is not, as a result of the bandpass. 

In our analysis, we construct an ensemble of statistics, each calculated by 
cross-correlating T = 60s worth of data for each instrument, which amounts to 
Ni = T/iStx) = 122880 points worth of LLO data and N 2 = T/(St 2 ) = 15000 points 
worth of ALLEGRO data. The LLO data are windowed and zero-padded out to a 
length Mi = 2N\ and discrete Fourier transformed to produce 

Nt-l Ni-1 

wSi[e}= wi\j]S!]j]e-W Ml = m[j}S 1 [j]e-^ e5 ^ 5t ^21) 
3=0 3=0 

where Sf = (Sh Ah)- 1 = (2T)- 1 = ^ Hz for T = 60s, and I = —Ni, ...,Ni~l. 
Meanwhile, the ALLEGRO data are windowed and zero-padded out to a length 
M 2 = 2N 2 and discrete Fourier transformed to produce 

JV3-1 

wS 2 [t] = W2[k}S^k}e- l2 ^ k / Ah e^ k5t2 

k=0 
N 2 -l 

= w 2 [k]S^[k]e- i2 < u f-^ k ^ (22) 

k=0 



Sf 



where 5f — (St 2 M 2 ) 1 is the same frequency resolution as before, and I — 

N 2 , ■ ■■ , ^jj +N 2 - 1. With these definitions, 

8tx, 2 ^S h2 [l]~Si, 2 il5f) (23) 

so if we construct a frequency-domain optimal filter Q(f), we can obtain an 
approximate analogue Y ~ Y c for the Y c defined in (fTTf) by constructing the statistic 

^raax 

Y= %(StiwSi\£\)*Q(l5f)(St2wS 2 [£\) (24) 

where (£ m i n Sf, £ max Sf) is the frequency range over which Q{f) has significant support. 
Because of the factor of P\{f)P 2 {f) in the denominator of l|14|) . this is limited to a 
subset of the frequency ranges associated with both the LLO and ALLEGRO data, 
as illustrated in Figure El 
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Figure 2. The frequency ranges covered by Fourier-transformed LLO and 
ALLEGRO data and the optimal filter used for analysis. The vertical scales 
are arbitrary and the data are meant to be illustrative only. If LLO data are 
downsampled to 2048 Hz, their discrete Fourier transform covers frequencies from 
-1024 Hz to 1024 Hz. The ALLEGRO data are heterodyned at 899 Hz and sampled 
at 250 Hz, which makes the frequencies represented range from 774 Hz to 1024 Hz. 
The optimal filter has non-negligible support only for 850 Hz < / < 950 Hz, so 
those are the only frequencies included in the cross-correlation. 



3.2. Time-Domain Equivalent 

In this section we consider a time-domain cross-correlation equivalent to the discrete 
frequency-domain approximation 124(1 and explicitly calculate the discrete-time 
equivalents of i|T2*|) and (|T3)l . 

One conceptually simple approach to cross-correlating these two data streams in 
the time domain would be to upsamplc both to the same sampling frequency. This 
should work as long as the ratio of the sampling rates is a rational number. We define 
ri and r 2 as the smallest integers such that 

St-t fi 

ot 2 r 2 

(in the present example, r\ = 125 and r 2 — 1024) and then define 

Jt=*^ = ^ (26) 
n r 2 

(In the present example, (St)^ 1 = 256000Hz.) These upsampled time series would 
then each contain 

N = Nin = N 2 r 2 (27) 

points. (In the present example, N = 15360000.) The method used to upsample 
the data does not interest us, since we won't actually perform the time-domain 
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upsampling. Instead, we idealize the result as the discrete sampling of the band- 
passed time series Si(t) and S^t): 

S[[J] =S 1 (t +j6t) J = 0,...,N (28) 
S% r [K] = S%(t +K8t) K = 0,...,N (29) 

If we also assume that + 1/2^2 < l/2St, there is no loss of information in going 
between the heterodyned and non-heterodyned data streams at the higher sampling 
rate. It is thus reasonable to think of the cross-correlation statistic as 

iV-l JV-l 

Y r = J2 6t l>2 St w l [ J l 5 i t J l *Q r i J - K ] w 2 [ R ] S 2 \ K \ (30) 

J=0 K=0 
JV-l JV-l 

= ^2 St E Stw ii J ] 5 iW Q r i J - K ) e l27r ^ K5t w r 2 [K] S% r [K] 

J=0 K=0 

(The factors of 5t are to facilitate comparison with the continuous-time idealization.) 
The optimal filter Q[J — K] depends only on the difference between the two indices, 
consistent with the assumption that we're considering stationary random processes, 
but the introduction of the two jV-point window functions§ wi^[J] (assumed to be 
real) allows us to control edge effects by smoothing out the onset and ending of the 
analyzed data, and the complex exponential accounts for the heterodyning. 

On the other hand, the analysis is not actually done with the upsampled time 
series Sl[J] and S2 r [K], but with the original S\[j] and S'^k]. So if we calculated a 
time-domain cross-correlation statistic, it would be 

JVi-l JV 2 -1 

Y = E E ^ 2 w ^ Q h ti> fc ] W2 w ^2 m (31) 

j=Q k=0 

Again, the bandpass tells us that the upsampled data streams don't have any higher- 
frequency content that's not present in the original ones. So we should get roughly 
the same cross-correlation statistic if we limit the sum in ()3U|I to J an integer multiple 
of n and K an integer multiple of r^, and then multiply by r\r2 to compensate for 
having taken fewer terms. This means Y r w Y with 

Q h [j, k] = Q r [r l3 - r 2 k] e l2 ^" k5t - (32) 

Since the sums over j and k in (|31|l both range from to N — 1, the argument 
°f Q[j ~ k] ranges from ~(N — 1) to N — 1, so a discrete Fourier transform (DFT) 
of Q will need to include at least 27V — 1 points. Since it is often more convenient to 
work with a 2iV-point DFT than a 2jV — 1-point one (e.g., if N is a power of two or a 
product of small primes), we will in general zero-pad Q[m] out to M > 2N — 1 points, 
with the "extra" values (i.e., those with N — 1 < m < M — 1) defined by 

0H-(° r , W - 1 <™<"-<"- 1 ) , (33) 

I Q[m -M] M-(N-l)<m<M v ' 

before defining the discrete Fourier transform 

M-l JV-l 

Qi £ ] = E e- t2 ™ e/M Q r [m] = e' t27rme/M Q r [m] . (34) 

m=0 m=-(JV-l) 

§ The analysis of windowing effects described here is a generalization of that described in 1171 and 
used in [7| for two data streams sampled at the same rate. 
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We can transform (|31[) into the frequency domain using the inverses of (@J) and (|34[l : 



S ia (t)= / dfe^f^S.M) 



M-l 



M ^ 



Q[i] 



(35) 
(36) 



£=0 



the result is 



£=0 W-oo 



([//-/i]T) -Si (A) 



4f2W2([//-/ a ]T)5 2 (/ a ) 

> 

where 5/ = (M(5t) _1 , /g = £5/, and the transformed window 



JV-l 



(37) 



(38) 



3=0 

is equivalent to an iVi or A^-point discrete Fourier transform, but not limited to 
integer arguments. Note that by construction Wi^(x) is periodic with period Ni^'- 
W 1>2 (x + N h2 ) = W 1 , 2 (x). 

3. 3. Mean and Variance of the Statistic 

We can get expressions for the expected mean and variance of (|37|l by applying © 
and 1|10|) . and noting that if we define Cij(f) by 

(si(frs j (f')) = S(f-f , )c ij (f) (39) 

so that cu(f) = \Pi{f) and c 12 (f) = |igw(/)7i2(/), tncn 

(Si(fyS j (f))=5(f-f)C ij (f) (40) 



where 



The mean is 



CM 



^(/) /2 l - 5 fe</<min( 5; l r ,^ + ^ 
otherwise 



(41) 



M-l 



=0 

oo 



d/ - /]T)* W 2 ([/, - f]T)C 12 {f) 



(42) 



As before, the restricted support of C%2(f) means that we can change the limits of 
the frequency integral from (—00,00) to (— jjt' TJt)' 

In this case, the periodicity of the windowing functions means that if they're 
sufficiently sharply peaked about zero argument, we must have (ft—f)T approximately 
equal to OmodTVi and 0modiV2. But since N is the lowest common multiple of JVi 
and N 2 , this is equivalent to the condition that (ft — f)T w OmodiV, which is true for 

at most one frequency in the range f£ — ^7 < / < mm fiSI' + ISt^j > an< ^ ^hat 
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frequency is always the positive fy. The windowing thus allows us to replace Cx%(f) 
with C\2{fl) hi (|42|l and obtain 



c min 

£ (Sti)(6h)Q[£]c 12 (f e ) 

dfW 1 ([f e -f]T)*W 2 ([f e -f}T) (43) 



M 

i—l • 

L l min 
rl/(2«) 



'-1/(2 r5t) 

where the limits on the sum over I are those such that 

f i-^<^<^{^ti + 2k) (44) 
The integral can be evaluated [using Q38J1] as 

4fWi([^-/]T)*^ 2 ([/^-/]T) 

-1/(2 at) 

Wi-lJV 2 -l A/(25t) 

= V V wi[j> 2 [fc] / e(f e&rUt-nUn-^st (45) 

! JVi-1^-1 1 ^ 

= 77 £ £ <W^lL?> 2 [fe] = 77—^1^2 

Note that the average of the product of the windows is taken only over the points for 
which both windows "coexist" given their different sampling rates: 

N/(rm) — l 

W1W2 = —~- ^ w 1 [nr 2 }w2[nri] (46) 

n=0 



This then tells us 



max 



H^wm^ ]T Sf(StQ[£\)-y(Jt)Pgw(ft) (47) 

c ^min 

where we have used again the definition 8f = l/(M St). 

As before, we can identify (|4Tjl (up to the factor of W1W2) as a discrete 
approximation to the usual continuous-time expression if we note that (|34[1 relates 
the discrete and continuous Fourier transforms according to ~ Q(fi). 

Note that if we design the filter in the frequency domain and chose Q[£] to be 
real, l|47|) tells us that the mean value of the statistic is real as long as any underlying 
correlation between the two data streams is time-symmetric. However, the band- 
passing means that s 2 (t) is complex, and therefore the statistic Y is as well. 

Since Y is complex, we consider the variance of x := KeY = Y+Y and 
y := ImY = separately. This ultimately comes down to calculating (Y 2 ) 

and (YY*), and specifically the contribution from the auto-correlations. 

In fact, the dominant contribution to (Y 2 ) vanishes because it contains 
(^(/O^Cf)}- Now, this is zero unless both / and /' lie in the range [f£ — ^Et^i ^2 + 
2gj^)i which in particular means both / and /' are positive. If it lies in that range, it 
is equal to 

(Mf)Mf)) = (M-f'TMf)) = Kf + f)^iP- (48) 
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but this vanishes unless / and /' have the opposite sign. Thus there is no combination 
of frequencies for which (S^C/')^!/)} is non-zero. 

This then means that the real and imaginary parts of Y both have a variance of 
1 rp M-i 2 

°i « °l = ^(Y*Y) ^-w\wl V \(StQ[£]) Pi(\M)P2(\h\)m 

where again the window average is conducted over the N/{r\r2) points where the 
windows "line up". 



4. Conclusions 



We have described the data analysis method used for cross-correlating in the 
frequency domain two data streams sampled at different rates, one of which is 
also heterodyned prior to digitization. We have illustrated heuristically how this 
approximates a continuous-time description, and also its equivalence to a hypothetical 
cross-correlation in the time domain using upsampled data streams. This frequency- 
domain technique allows us to efficiently cross-correlate data from different detectors 
with different sampling parameters, and is being applied to LIGO and ALLEGRO 
data taken during LIGO's second science run. 
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